Response Time-dependent Density Functional Theory Research Articles

A mathematical analysis of the adiabatic Dyson equation from time-dependent density functional theory

In this article, we analyse the Dyson equation for the density–density response function (DDRF) that plays a central role in linear response time-dependent density functional theory (LR-TDDFT). First, we present a functional analytic setting that allows for a unified treatment of the Dyson equation with general adiabatic approximations for discrete (finite and infinite) and continuum systems. In this setting, we derive a representation formula for the solution of the Dyson equation in terms of an operator version of the Casida matrix. While the Casida matrix is well-known in the physics literature, its general formulation as an (unbounded) operator in the N-body wavefunction space appears to be new. Moreover, we derive several consequences of the solution formula obtained here; in particular, we discuss the stability of the solution and characterise the maximal meromorphic extension of its Fourier transform. We then show that for adiabatic approximations satisfying a suitable compactness condition, the maximal domains of meromorphic continuation of the initial DDRF and the solution of the Dyson equation are the same. The results derived here apply to widely used adiabatic approximations such as (but not limited to) the random phase approximation and the adiabatic local density approximation. In particular, these results show that neither of these approximations can shift the ionisation threshold of the Kohn–Sham system.

Prediction of Strong Solvatochromism in a Molecular Photocatalyst.

Based on quantum chemical calculations, we predict strong solvatochromism in a light-driven molecular photocatalyst for hydrogen generation, that is we show that the electronic and optical properties of the photocatalyst strongly depend on the solvent it is dissolved in. Our calculations in particular indicate a solvent-dependent relocation of the highest occupied molecular orbital (HOMO). Ground-state density functional theory and linear response time-dependent density functional theory calculations were applied in order to investigate the influence of implicit solvents on the structural, electronic and optical properties of a molecular photocatalyst. Only at high dielectric constants of the solvent, is the HOMO located at the metal center of the photosensitizer, whereas at low dielectric constants the HOMO is centered at the metal atom of the catalytically active complex. We elucidate the electronic origins of this strong solvatochromic effect and sketch the consequences of these insights for the use of photocatalysts in different environments.

Designing Efficient Metal-Free Dye-Sensitized Solar Cells: A Detailed Computational Study.

The modulation of molecular characteristics in metal-free organic dyes holds significant importance in dye-sensitized solar cells (DSSCs). The D-π-A molecular design, based on the furan moiety (π) in the conjugated spacer between the arylamine (D) and the 2-cyanoacrylic acid (A), was developed and theoretically evaluated for its potential application in DSSCs. Utilizing linear response time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional, different donor and acceptor groups were characterized in terms of the electronic absorption properties of these dyes. All the studied dye sensitizers demonstrate the ability to inject electrons into the semiconductor's conduction band (TiO2) and undergo regeneration through the redox potential triiodide/iodide (I3-/I-) electrode. TDDFT results indicate that the dyes with CSSH anchoring groups exhibit improved optoelectronic properties compared to other dyes. Further, the photophysical properties of all dyes absorbed on a Ti(OH)4 model were explored and reported. The observed results indicate that bidentate chemisorption occurs between dyes and TiO4H5. Furthermore, the HOMO-LUMO energy gaps for almost all dye complexes are significantly smaller than those of the free dyes. This decrease of the HOMO-LUMO energy gaps in the dye complexes facilitates electron excitation, and thus more photons can be adsorbed, guaranteeing larger values of efficiency and short-circuit current density.

Density Functional Treatment of Photoionization of Sodium Clusters: Effects of Cluster Size and Exchange–Correlation Framework

The ground state and photoionization properties of Nax (x = 20, 40, and 92) clusters are investigated using a method based on density functional theory (DFT) in a spherical jellium frame. Two different exchange–correlation treatments with the Gunnarsson–Lundqvist parametrization are used: (i) the electron self-interaction correction (SIC) scheme and (ii) the van Leeuwen–Baerends (LB94) scheme based on the gradient of the electron density. The shapes of the mean-field potentials and bound state properties, obtained in the two schemes, qualitatively agree, but differ in the details. The effect of the schemes on the photoionization dynamics, calculated in linear response time-dependent DFT is compared, in which the broader features are found to be universal. The general similarity of the results in SIC and LB94 demonstrates the reliability of DFT treatments. The study further elucidates the evolution of the ground state and ionization description as a function of the cluster size.

Core-Projected Hybrids Fix Systematic Errors in Time-Dependent Density Functional Theory Predicted Core-Electron Excitations.

Linear response time-dependent density functional theory (TDDFT) is widely applied to valence, Rydberg, and charge-transfer excitations but, in its current form, makes large errors for core-electron excitations. This work demonstrates that the admixture of nonlocal exact exchange in atomic core regions significantly improves TDDFT-predicted core excitations. Exact exchange admixture is accomplished using projected hybrid density functional theory [ J. Chem. Theory Comput. 2023, 19, 837-847]. Scalar relativistic TDDFT calculations using core-projected B3LYP accurately model core excitations of second-period elements C-F and third-period elements Si-Cl, without sacrificing performance for the relative shifts of core excitation energies. Predicted K-edge X-ray near absorption edge structure (XANES) of a series of sulfur standards highlight the value of this approach. Core-projected hybrids appear to be a practical solution to TDDFT's limitations for core excitations, in the way that long-range-corrected hybrids are a practical solution to TDDFT's limitations for Rydberg and charge-transfer excitations.

Core Excitations of Uranyl in Cs2UO2Cl4 from Relativistic Embedded Damped Response Time-Dependent Density Functional Theory Calculations.

X-ray spectroscopies, by their high selectivity and sensitivity to the chemical environment around the atoms probed, provide significant insights into the electronic structures of molecules and materials. Interpreting experimental results requires reliable theoretical models, accounting for environmental, relativistic, electron correlation, and orbital relaxation effects in a balanced manner. In this work, we present a protocol for the simulation of core excited spectra with damped response time-dependent density functional theory based on the Dirac-Coulomb Hamiltonian (4c-DR-TD-DFT), in which environmental effects are accounted for through the frozen density embedding (FDE) method. We showcase this approach for the uranium M4- and L3-edges and oxygen K-edge of the uranyl tetrachloride (UO2Cl42-) unit as found in a host Cs2UO2Cl4 crystal. We have found that the 4c-DR-TD-DFT simulations yield excitation spectra that very closely match the experiment for the uranium M4-edge and the oxygen K-edge, with good agreement for the broad experimental spectra for the L3-edge. By decomposing the complex polarizability in terms of its components, we have been able to correlate our results with angle-resolved spectra. We have observed that for all edges, but in particular the uranium M4-edge, an embedded model in which the chloride ligands are replaced by an embedding potential reproduces rather well the spectral profile obtained for UO2Cl42-. Our results underscore the importance of the equatorial ligands to simulating core spectra at both uranium and oxygen edges.

Analytical excited state gradients for time-dependent density functional theory plus tight binding (TDDFT + TB).

Understanding photoluminescent mechanisms has become essential for photocatalytic, biological, and electronic applications. Unfortunately, analyzing excited state potential energy surfaces (PESs) in large systems is computationally expensive, and hence limited with electronic structure methods such as time-dependent density functional theory (TDDFT). Inspired by the sTDDFT and sTDA methods, time-dependent density functional theory plus tight binding (TDDFT + TB) has been shown to reproduce linear response TDDFT results much faster than TDDFT, particularly in large nanoparticles. For photochemical processes, however, methods must go beyond the calculation of excitation energies. Herein, this work outlines an analytical approach to obtain the derivative of the vertical excitation energy in TDDFT + TB for more efficient excited state PES exploration. The gradient derivation is based on the Z vector method, which utilizes an auxiliary Lagrangian to characterize the excitation energy. The gradient is obtained when the derivatives of the Fock matrix, the coupling matrix, and the overlap matrix are all plugged into the auxiliary Lagrangian, and the Lagrange multipliers are solved. This article outlines the derivation of the analytical gradient, discusses the implementation in Amsterdam Modeling Suite, and provides proof of concept by analyzing the emission energy and optimized excited state geometry calculated by TDDFT and TDDFT + TB for small organic molecules and noble metal nanoclusters.

Calculations of Excited Electronic States by Converging on Saddle Points Using Generalized Mode Following.

Calculations of excited electronic states are carried out by finding saddle points on the surface describing the variation of the energy of the system as a function of the electronic degrees of freedom. This approach has several advantages over commonly used methods especially in the context of density functional calculations, as collapse to the ground state is avoided, and yet, the orbitals are variationally optimized for the excited state. Such a state-specific optimization makes it possible to describe excitations with large charge transfer, where calculations based on ground state orbitals, such as linear response time-dependent density functional theory, can be problematic. A generalized mode following method is presented where an nth-order saddle point is found by inverting the components of the gradient in the direction of the eigenvectors of the n lowest eigenvalues of the electronic Hessian matrix. This approach has the distinct advantage of following a chosen excited state by its saddle point order through molecular configurations where the symmetry of the single determinant wave function is broken, thereby making it possible to calculate potential energy curves even at avoided crossings, as demonstrated here in calculations of the ethylene and dihydrogen molecules. Results of calculations are, furthermore, presented for charge transfer excitations in nitrobenzene and N-phenylpyrrole, corresponding to fourth- and sixth-order saddle points, respectively, where an approximate initial estimate of the saddle point order could be found by energy minimization with excited electron and hole orbitals frozen. Finally, calculations of a diplatinum-silver complex are presented, illustrating the applicability of the method to larger molecules.

Curing the Divergence in Time-Dependent Density Functional Quadratic Response Theory

The adiabatic approximation in time-dependent density functional theory is known to give an incorrect pole structure in the quadratic response function, leading to unphysical divergences in excited state-to-state transition probabilities and hyperpolarizabilties. We find the form of the exact quadratic response kernel and derive a practical and accurate approximation that cures the divergence. We demonstrate our results on excited state-to-state transition probabilities of a model system and of the LiH molecule.

Theoretical Investigations on the Plasmon-Mediated Dissociation of Small Molecules in the Presence of Silver Atomic Wires.

Plasmonic nanoparticles can promote bond activation in adsorbed molecules under relatively benign conditions via excitation of the nanoparticle's plasmon resonance. As the plasmon resonance often falls within the visible light region, plasmonic nanomaterials are a promising class of catalysts. However, the exact mechanisms through which plasmonic nanoparticles activate the bonds of nearby molecules are still unclear. Herein, we evaluate Ag8-X2 (X = N, H) model systems via real-time time-dependent density functional theory (RT-TDDFT), linear response time-dependent density functional theory (LR-TDDFT), and Ehrenfest dynamics in order to better understand the bond activation processes of N2 and H2 facilitated by the presence of the atomic silver wire under excitation at the plasmon resonance energies. We find that dissociation is possible for both small molecules at high electric field strength. Activation of each adsorbate is symmetry- and electric field-dependent, and H2 activates at lower electric field strengths than N2. This work serves as a step toward understanding the complex time-dependent electron and electron-nuclear dynamics between plasmonic nanowires and adsorbed small molecules.

State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives.

Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au25(SR)18 - is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.

Exact Two-Component TDDFT with Simple Two-Electron Picture-Change Corrections: X-ray Absorption Spectra Near L- and M-Edges of Four-Component Quality at Two-Component Cost.

X-ray absorption spectroscopy (XAS) has gained popularity in recent years as it probes matter with high spatial and elemental sensitivities. However, the theoretical modeling of XAS is a challenging task since XAS spectra feature a fine structure due to scalar (SC) and spin-orbit (SO) relativistic effects, in particular near L and M absorption edges. While full four-component (4c) calculations of XAS are nowadays feasible, there is still interest in developing approximate relativistic methods that enable XAS calculations at the two-component (2c) level while maintaining the accuracy of the parent 4c approach. In this article we present theoretical and numerical insights into two simple yet accurate 2c approaches based on an (extended) atomic mean-field exact two-component Hamiltonian framework, (e)amfX2C, for the calculation of XAS using linear eigenvalue and damped response time-dependent density functional theory (TDDFT). In contrast to the commonly used one-electron X2C (1eX2C) Hamiltonian, both amfX2C and eamfX2C account for the SC and SO two-electron and exchange-correlation picture-change (PC) effects that arise from the X2C transformation. As we demonstrate on L- and M-edge XAS spectra of transition metal and actinide compounds, the absence of PC corrections in the 1eX2C approximation results in a substantial overestimation of SO splittings, whereas (e)amfX2C Hamiltonians reproduce all essential spectral features such as shape, position, and SO splitting of the 4c references in excellent agreement, while offering significant computational savings. Therefore, the (e)amfX2C PC correction models presented here constitute reliable relativistic 2c quantum-chemical approaches for modeling XAS.

Automatic purpose-driven basis set truncation for time-dependent Hartree–Fock and density-functional theory

Real-time time-dependent density-functional theory (RT-TDDFT) and linear response time-dependent density-functional theory (LR-TDDFT) are two important approaches to simulate electronic spectra. However, the basis sets used in such calculations are usually the ones designed mainly for electronic ground state calculations. In this work, we propose a systematic and robust scheme to truncate the atomic orbital (AO) basis set employed in TDDFT and TD Hartree–Fock (TDHF) calculations. The truncated bases are tested for both LR- and RT-TDDFT as well as RT-TDHF approaches, and provide an acceleration up to an order of magnitude while the shifts of excitation energies of interest are generally within 0.2 eV. The procedure only requires one extra RT calculation with 1% of the total propagation time and a simple modification on basis set file, which allows an instant application in any quantum chemistry package supporting RT-/LR-TDDFT calculations. Aside from the reduced computational effort, this approach also offers valuable insight into the effect of different basis functions on computed electronic excitations and further ideas on the design of basis sets for special purposes.

Structural and spectroscopic properties of voriconazole and fluconazole – Experimental and theoretical studies

Abstract The article compares the experimental Fourier transform-infrared, Ultraviolet-visible (UV-vis), and proton nuclear magnetic resonance (1H NMR) spectra of voriconazole and fluconazole with the density functional theory (DFT) calculations using five different functionals. The results were compared with previously reported data related to its analogue – posaconazole. The analysis of calculated infrared (IR) spectra with the use of PBE1PBE (voriconazole) or APF (fluconazole) functionals shows good accordance with the experimental IR spectrum. The best compatibility between the experimental and theoretical UV spectra was observed with the use of PBE1PBE or B3LYP functionals for voriconazole or fluconazole, respectively. The reason for the difference in the UV-vis spectra of voriconazole and fluconazole was discussed based on linear response time-dependent DFT and natural bond orbital methods. The calculated 1H NMR spectrum shows that the DFT formalism, particularly the M06L or B3LYP functionals, gives an accurate description of the voriconazole and fluconazole chemical shifts.

Structural and Spectroscopic Properties of Isoconazole and Bifonazole-Experimental and Theoretical Studies.

The paper compares the experimental FT-IR, UV-vis, and 1H NMR spectra of isoconazole and bifonazole with the density functional theory (DFT) calculations using different functionals. The results were compared with previously reported data related to their analogue, posaconazole. The analysis of calculated IR spectra with use of CAM-B3LYP (isoconazole) or B3LYP (bifonazole) functionals shows good accordance with the experimental IR spectrum. The best compatibility between the experimental and theoretical UV spectra was observed with the use of B3LYP or wB97XD functionals for isoconazole or bifonazole, respectively. The reason for the difference in the UV-vis spectra of isoconazole and bifonazole was discussed based on linear response time-dependent DFT and natural bond orbital methods. The calculated 1H NMR spectrum shows that the DFT formalism, particularly the B3LYP functional, give an accurate description of the isoconazole and bifonazole chemical shifts.

Investigation of plasmon relaxation mechanisms using nonadiabatic molecular dynamics.

Hot carriers generated from the decay of plasmon excitation can be harvested to drive a wide range of physical or chemical processes. However, their generation efficiency is limited by the concomitant phonon-induced relaxation processes by which the energy in excited carriers is transformed into heat. However, simulations of dynamics of nanoscale clusters are challenging due to the computational complexity involved. Here, we adopt our newly developed Trajectory Surface Hopping (TSH) nonadiabatic molecular dynamics algorithm to simulate plasmon relaxation in Au20 clusters, taking the atomistic details into account. The electronic properties are treated within the Linear Response Time-Dependent Tight-binding Density Functional Theory (LR-TDDFTB) framework. The relaxation of plasmon due to coupling to phonon modes in Au20 beyond the Born-Oppenheimer approximation is described by the TSH algorithm. The numerically efficient LR-TDDFTB method allows us to address a dense manifold of excited states to ensure the inclusion of plasmon excitation. Starting from the photoexcited plasmon states in Au20 cluster, we find that the time constant for relaxation from plasmon excited states to the lowest excited states is about 2.7ps, mainly resulting from a stepwise decay process caused by low-frequency phonons of the Au20 cluster. Furthermore, our simulations show that the lifetime of the phonon-induced plasmon dephasing process is ∼10.4fs and that such a swift process can be attributed to the strong nonadiabatic effect in small clusters. Our simulations demonstrate a detailed description of the dynamic processes in nanoclusters, including plasmon excitation, hot carrier generation from plasmon excitation dephasing, and the subsequent phonon-induced relaxation process.

Effect of indium vacancy point defect on the structural, electronic and optical properties of double perovskite halide Cs2InSbCl6

Efforts to replace the lead-based perovskite have gain significant improvements in the last decade. Lead free double halide perovskite are being considered to replace the toxic lead-based perovskite materials for solar cell applications. Density functional theory (DFT) and linear response time-dependent density functional theory (TDDFT) are used in this study to simulate and investigate the effect of vacancy in the Indium (In) site of leadfree double halide perovskite A2BB′X6 (A = Cs; B = In, B′ = Sb; X = Cl) on the structural, electronic, and optical properties for possible solar cell application. On the structural properties, the total bond length of the material Cs2InSbCl6 was found to be 30.3848 Å and the creation of the indium vacancy raises the value to 30.594 Å. The bulk modulus of Cs2InSbCl6 was calculated to be 28.44 Gpa and 26.22 Gpa for Cs2In-xSbCl6. The calculated band structure for Cs2InSbCl6 reveals semiconducting behavior with a direct energy band gap of 0.99 eV along 𝛤 point symmetry while the created indium vacancy band structure reveals metallic behavior with an overlap of electronic states along the Γ - W – Γ symmetry. The result of the optical properties calculations shows that the material Cs2InSbCl6 has higher absorption coefficient, low refractive index, low reflectivity and higher conductivity when compared to the created indium vacancy material.

Orbital-free methods for plasmonics: Linear response.

Plasmonic systems, such as metal nanoparticles, are widely used in different areas of application, going from biology to photovoltaics. The modeling of the optical response of such systems is of fundamental importance to analyze their behavior and to design new systems with required properties. When the characteristic sizes/distances reach a few nanometers, nonlocal and spill-out effects become relevant and conventional classical electrodynamics models are no more appropriate. Methods based on the Time-Dependent Density Functional Theory (TD-DFT) represent the current reference for the description of quantum effects. However, TD-DFT is based on knowledge of all occupied orbitals, whose calculation is computationally prohibitive to model large plasmonic systems of interest for applications. On the other hand, methods based on the orbital-free (OF) formulation of TD-DFT can scale linearly with the system size. In this Review, OF methods ranging from semiclassical models to the Quantum Hydrodynamic Theory will be derived from the linear response TD-DFT, so that the key approximations and properties of each method can be clearly highlighted. The accuracy of the various approximations will then be validated for the linear optical properties of jellium nanoparticles, the most relevant model system in plasmonics. OF methods can describe the collective excitations in plasmonic systems with great accuracy and without system-tuned parameters. The accuracy of these methods depends only on the accuracy of the (universal) kinetic energy functional of the ground-state electronic density. Current approximations and future development directions will also be indicated.

Accurate Calculation of Excited-State Absorption for Small-to-Medium-Sized Conjugated Oligomers: Multiconfigurational Treatment vs Quadratic Response TD-DFT.

Excited-state absorption (ESA) spectra of π-conjugated compounds are frequently calculated by (quadratic response) time-dependent density functional theory, (QR) TD-DFT, often giving a reasonable representation of the experimental results despite the (known) incomplete electronic description. To investigate whether this is inherent to the method, we calculate here the ESA spectra of small-to-medium-sized oligophenylenevinylenes (nPV) and oligothiophenes (nT) using QR TD-DFT as well as CASPT2 based on CASSCF geometries. CASPT2 gives indeed a reliable, theoretically correct description of the ESA features for all compounds; the computational effort can be reduced without significant loss of accuracy using TD-DFT geometries. QR TD-DFT, based on BHandHLYP and CAM-/B3LYP functionals, fails on short nTs but provides a reasonable description for spectral positions of nPVs and long nTs. The failure on short nTs is, however, only partly due to the incomplete configuration description but, in particular, related to an improper MO description, resulting in an asymmetric energy spacing of the occupied vs unoccupied MOs in the DFT scheme. Longer nTs, on the other side, adapt approximately the MO scheme for alternant hydrocarbons just like in nPVs, while contributions by two triplet excitations combined to a singlet (which inhibits an accurate treatment of polyenes with standard TD-DFT) do not play a relevant role in the current case. For such "well-behaved" systems, a reasonable representation of ESA spectra is found at the QR TD-DFT level due to the rather small energy shifts when including higher-order excitations.

High-frequency limit of spectroscopy.

We consider an arbitrary quantum mechanical system, initially in its ground-state, exposed to a time-dependent electromagnetic pulse with a carrier frequency ω0 and a slowly varying envelope of finite duration. By working out a solution to the time-dependent Schrödinger equation in the high-ω0 limit, we find that, to the leading order in ω0 -1, a perfect self-cancellation of the system's linear response occurs as the pulse switches off. Surprisingly, the system's observables are, nonetheless, describable in terms of a combination of its linear density response function and nonlinear functions of the electric field. An analysis of a jellium slab and jellium sphere models reveals a very high surface sensitivity of the considered setup, producing a richer excitation spectrum than accessible within the conventional linear response regime. On this basis, we propose a new spectroscopic technique, which we provisionally name the Nonlinear High-Frequency Pulsed Spectroscopy (NLHFPS). Combining the advantages of the extraordinary surface sensitivity, the absence of constraints by the traditional dipole selection rules, and the clarity of theoretical interpretation utilizing the linear response time-dependent density functional theory, NLHFPS has a potential to evolve into a powerful characterization method for nanoscience and nanotechnology.